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ABSTRACT 

We analyze star formation (SF) as a function of stellar mass (M*) and redshift z in the All Wave- 
length Extended Groth Strip International Survey (AEGIS), for star-forming field galaxies with 
M* > 1O 1O M out to z = 1.1. The data indicate that the high specific SF rates (SFR) of many 
less massive galaxies do not represent late, irregular or recurrent, starbursts in evolved galaxies. They 
rather seem to reflect the onset (initial burst) of the dominant SF episode of galaxies, after which SF 
gradually declines on Gyr timescales to z = and forms the bulk of a galaxy's M*. With decreasing 
mass, this onset of major SF shifts to decreasing z for an increasing fraction of galaxies (staged galaxy 
formation). This process may be an important component of the "downsizing" phenomenon. We 
find that the predominantly gradual decline of SFR described by Noeske et al. can be reproduced by 
exponential SF histories (t models) , if less massive galaxies have systematically longer e- folding times 
r, and a later onset of SF (zi). Our model can provide a first parametrization of SFR as a function 
of M+ and z, and quantify mass-dependences of r and Zf, from direct observations of M+ and SFR up 
to z > 1. The observed evolution of SF in galaxies can plausibly reflect the dominance of gradual gas 
exhaustion. The data are also consistent with the history of cosmological accretion onto Dark Matter 
halos. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: 
starburst 



1. INTRODUCTION 

In an accompanying Letter (Noeske et al. 2007, here- 
after Paper I), we have studied star formation (SF) rates 
(SFR) as a function of stellar mass (M+) and z, for field 
galaxies in the All Wavelength Extended Groth Strip 
International Survey (AEGIS) out to z = 1.1. Star- 
forming galaxies form a defined "Main Sequence" (MS) 
with a limited range of SFR at a given M* and z. This 
smooth sequence suggests that the same set of few phys- 
ical processes governs SF in galaxies, unless quenching 
of their SF occurs. The evolution of SF along the MS 
appears dominated by a gradual decline of SFR in indi- 
vidual galaxies since z ~ 1 , not by an evolving frequency 
or amplitude of starbursts. The dominant process that 
governs SF since z ~ 1 is hence likely a gradual one, an 
obvious candidate being gas exhaustion. 

SF histories (SFHs) are known to depend on galaxy 
mass and morphological type, both from studies of local 
galaxies (e.g. Tinsley 1968, Searle et al. 1973, Sandage 
1986, Heavens et al. 2004), and from distant galaxy sur- 
veys (see references in §1 of Paper I). The common pic- 
ture is that massive galaxies formed the bulk of their 
stars early and on shorter timescales, while numerous 
less massive galaxies evolve on longer timescales, a phe- 
nomenon generally linked to the "downsizing" reported 
by Cowie et al. (1996). 
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In this Letter, we show that a simple model of gas ex- 
haustion with mass-dependent parameters can reproduce 
and parametrize the observed SFR as a function of M* 
and z. Gas exhaustion may thus be responsible for the 
gradual decline of SFR that dominates SFHs since z ~ 1 
along the MS, i.e. in star-forming field galaxies. Follow- 
ing previous authors, we consider specific SFR (SSFR), 
i.e. SFR/M*, a simple but powerful indicator of galaxy 
SFHs (e.g. Kennicutt et al. 2005). We argue that the 
onset of major SF occurs systematically later in less mas- 
sive galaxies. 

2. DATA SET 

As in Paper I, we take advantage of the sensitivity and 
panchromatic nature of AEGIS; combined SFRs from 
deep Spitzer MIPS 24/im images and DEEP2 spectra re- 
cover obscured SF in IR-luminous galaxies, and achieve 
a large dynamic range in SFR by including galaxies not 
detected at 24/iin. For a description of the data, SFR 
tracers and JW* measurements, see Section 2 of Paper I. 
We consider all galaxies with robust SFR tracers a MS 
galaxy - either 24/zm detected, or blue sequence galax- 
ies with S/N > 2 emission lines (Ha, H/3, or [OII]3727), 
thereby excluding red LINER/ AGN candidates (see Pa- 
per I). As shown in Paper I, this selection likely misses at 
most < 10(20)% of the normally star-forming MS galax- 
ies at z < (>)0.7, likely less. 

We tested the effects of using different combinations 
of SFR and M* measures, including GALEX UV based 
SFRs, and M* values from the color-M/L relation of Bell 
et al (2003). All qualitative results of this work are robust 
against the choice of SFR tracer or M+ estimate, yet 
quantitative results will vary (see §4.1). 

3. PARAMETRIZATION THROUGH r MODELS 
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Interpreting the SSFR vs {M+, z) diagrams in terms 
of mass-dependent SFH is not straightforward, as M* 
grows with time for SF galaxies. Here we present the 
use of a simple exponential model SFH (r models, Eq. 
1) with mass-dependent parameters to quantify mass 
dependences of SF timescales, and to account for M* 
growth. Previous authors have successfully employed r 
models with different e-folding times r to reproduce the 
spectrophotomctric and chemical evolution of different 
Hubble types and masses (e.g. Tinsley 1968, Searle et 
al. 1973, Koo et al. 1993, Bicker et al. 2004, Savaglio 
et al. 2005, Weiner et al. 2006). The apparent dom- 
inance of smoothly declining SFR in individual galax- 
ies (Paper I) supports the use of r models, which are 
a one-zone approach to describe SF through continuous 
gas exhaustion. We adopt simple closed-box conditions 
where galaxies have a baryonic mass Mj that is initially 
gaseous, later the sum of gas (M g ) and stellar mass M*. 
For instantaneous recycling, with a recycled gas fraction 
R = 0.5 (Kroupa IMF, Bell et al. 2005), and a SF effi- 



ciency e such that the SFR \& = eM g , one obtains: 



*(M b ,z) = *(z/)exp 
T = t(z)-t(z f ),r = 
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1 



e(l-R) 



(1) 
(2) 



where Zf is the "formation redshift" where SF begins and 
t(z) is the cosmic time at redshift z The initial SFR at 
a given r is then = eM b = [t(1 - R)]~ 1 M b . We 

parametrize the mass dependence of r as a power law of 
the baryonic mass of the galaxy M b : 



r(M b 



(3) 



Fig. 1 (left column) shows examples of Eq. 1 in the 
SSFR-Af* plane, compared to the median SSFR of the 
MS, for different zj, c a and a. 

3.1. Staged r models 

Fig. 1 (left column) shows that models with mass- 
dependent r can crudely reproduce the median MS of SF 
galaxies and its redshift evolution for galaxies with M* > 
10 10 Af Q out to z ~ 1, if formation redshifts Zf ~ 2 are 
adopted for all galaxies. However, the models fall short of 
reproducing the high SSFR of less massive galaxies. The 
model SSFR remain systematically too low, unless we 
adopt a very low Zf < 1, unphysical for massive galaxies. 
The reason is the monotonic decline of the SFR of t 
models. Their present-to-past average SFR (Kennicutt 
et al. 2005), 



b(t) 



< * > T ~~ M* 1 — R' 



(4) 



is always < 1. The limit for r = oo is b = 1, which cor- 
responds to the constant SFR that would have formed a 
galaxy's M+ since Zf. Empirically, the behavior of the 
MS suggests declining SFHs (b < 1), which causes a con- 
flict between the high SSFR for low mass galaxies and 
the assumption that all galaxies started forming stars 
at high ZfQ> 3). Adopting such high Zf, an early start 
of SF, implies low past-average SFR: 6=1 then corre- 
sponds to low SSFR, reflected in the low upper limits 
to the SSFR of r models. For these Zf, the high SSFR 



of many less massive galaxies would imply b > 1, a cur- 
rent SFR above the past-average level, i.e. an episode of 
enhanced SF (Fig. 1, left; Fig. 2). 

The increase of the highest SSFR towards less massive 
galaxies can be reproduced by allowing the onset of SF 
to be delayed to lower zj in less massive galaxies (see the 
models for different Zf in Fig. 1, left). We parametrize 
zj as a function of mass, similarly to t, an approach we 
refer to as staged t models. 

l + z f (M b )=cpM b 9 (5) 

This model interprets high SSFR as the early epoch of 
smooth SFHs with lower Zf, rather than late episodes of 
enhanced SF. It is physically motivated, not an attempt 
to force an oversimplified model to fit complex SFHs (see 
§4.2). Staged t models (right column in Fig. 1) provide a 
better description of the median of the MS than r models 
with fixed Zf in the M + range where the sample is com- 
plete, and appear to describe the data also towards lower 
M*. Staged models that consider a moderate range of Zf 
and t at a given mass (dashed red curves) also reproduce 
the upper envelope of the MS, which is complete at all 
observed M* and z, and the lower envelope and apparent 
broadening of the MS towards lower masses. Models with 
a range of Zf but no trend with mass would merely intro- 
duce a scatter and an offset in the asymptotic SSFR at 
low masses, but would not change the slope of SSFR(Af^) 
to first order (see the models for different z/ in Fig. 1, 
left column). 
The staged r model in Fig. 1 (right) is given by 
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( 6 ) 

We calculated % 2 to scan the model parameter space, but 
Eq. [5] is hand-adjusted to reproduce the median and up- 
per and lower limits of the data. Results of a simple % 2 
minimization to the median would be misleading: best 
values for all four parameters in Eqs. (3) and (5) de- 
pend considerably on systematics of, e.g., SFR and M* 
estimates, and the IMF; also, the r(M b ) dependence is 
mainly constrained by massive galaxies (Fig. 1) where 
our number statistics are poor, and the Zf(M b ) relation 
at low masses, where data are incomplete. An evaluation 
of the relevant uncertainties must incorporate scatter in 
r(M b ) and Zf(M b ), or a scatter about smooth SF histo- 
ries, and is postponed to a forthcoming paper. 

4. DISCUSSION 

4.1. t models, gradual decline of SF 

By direct measurement of SFR and M* over a large 
range in mass and z, we confirm that the commonly 
adopted exponential model SFHs can reproduce the av- 
erage SFH of MS galaxies. This model can quantify the 
mass dependence of the associated SF timescales r and 
of the Zf. The mass dependences of r and Zf, and M+ 
growth through SF, conspire to reproduce the decline of 
SFR that is similar over a wide M* range (Paper I; Zheng 
et al. 2007). 

Notably, r models are a simple approximation of SF 
that declines due to gradual gas exhaustion. Their abil- 
ity to reproduce the evolution of the MS of SF galax- 
ies, along with the limited range of SFR on the MS, im- 
plies that gradual gas exhaustion with mass-dependent 
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Fig. 1. — Specific SFR [yr~ 1 ] vs M* for 3658 star-forming ("Main Sequence", Paper I) AEGIS galaxies. Full blue dots: SFR from Spitzer 
MIPS 24/im and DEEP2 emission lines (Paper I); open blue dots: blue galaxies without 24/im detection, SFR from extinction-corrected 
emission lines. Galaxies with no reliable signs of SF, inch red LINER/AGN candidates (Paper I), are not shown. Black dots and error bars: 
median and sample standard deviation of log(SSFR) of the main sequence galaxies, in the M* range where the sample is > 95% complete. 
The black dot-dashed line repeats the green (left) and red (right) models in the lowest z bin. left: r models with fixed formation redshift zj 
and mass-dependent r (colored curves). Massive galaxies can be reproduced assuming high z/(> 2), less massive galaxies require Zf < 2, 
unphysical for massive galaxies, right: Staged r models (red), where both r and Zf are mass-dependent. Red dashed lines show the effect 
of varying Zf and r at a given M^ aT . The delayed onset of SF (lower Zf) in a fraction of less massive galaxies accounts for the increase of 
specific SFR at low JW* without requiring a large fraction of galaxies to simultaneously undergo starbursts. 



timescales is a plausible driver of the dominant evolution 
of SF in galaxies > 1O 1O M0 since z ~ 1. 

We chose a closed box model which is sufficient to 
reproduce the coevolution of SF and M*. Linking the 
model M* to, e.g. Dark Matter halo masses should in- 
volve the observed relation between both values at a 
given z, to account for gas accretion and removal in 
galaxies, which are both not well understood. The r 
models' similarity to the data does not imply a closed 
box scenario where gas is merely turned into stars. Addi- 
tional processes that gradually deplete cold gas — heat- 
ing or loss — and scale roughly with SFR would also pro- 
duce SFHs that resemble exponentials. These processes 
include feedback from SF, and conceivably AGN, given 
the likely co-evolution of stellar bulges and Black Holes 
(e.g. Granato et al. 2004). Short r obtained for massive 
galaxies may largely reflect such gas loss processes rather 
than very efficient SF. 

We have considered the depletion of an existing gas 
reservoir, but the decline of SFR at a given mass is 
also compatible with the cosmological decline in accre- 
tion onto dark halos. This can be approximated for halo 
masses near 1O 12 M by M h /M h w 0.04 Gyr" 1 (1 + z) 2 25 



(Birnboim et al. 2007), giving a factor of ~ 5 between 
z ~ 1 and z ~ 0, similar to the observed decline. The 
mean virial accretion in this mass range is predicted to 

vary as M h oc M^ 15 . If we adopt M+fMh oc V 2 oc M^ 3 
(Dekel & Woo 2003 for SN feedback), and, naively, 
M* oc Mh, we obtain M* oc Aff 69 , compatible with the 
observed mass dependence (Paper I). 

4.2. Staged galaxy formation 

Fig. 1 shows the previously described increase of the 
highest SSFR towards lower M*. High SSFR have been 
interpreted as episodes of SFR above the past average 
level (b > 1), based on the assumption that all galaxies 
had a relatively high z/ (see §3) (e.g., Bell et al. 2005, 
Juneau et al. 2005, Feulner et al. 2005). 

Our data indicate that high SSFR often do not repre- 
sent such irregular or periodic episodes late in the SFH 
of a galaxy. In Fig. 2, we show SSFR vs M* as in 
Fig. 1. We express SSFR in terms of the doubling times, 
t d = M+/[(l - i?)*] = [(1 - R) SSFR]- 1 , within which 
the current SFR would produce the observed M*. Small 
td corresponds t high SSFR; for a declining SFH, a galaxy 
can be at most as old as td- Consider the galaxies at 
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Fig. 2. — Stellar mass doubling times (see §4.2) as a function 
of Mi,. Symbols are as in Figure 1. Horizontal lines give the time 
between Zf and the center of the z bin shown, for Zf of 1 and 5. 
Galaxies above these lines would have shorter than their adopted 
age, i.e. SFR above their past-average level (b > 1). Galaxies on 
the lines form stars at their past-average rate (b = 1). 

0.75 < z < 1 and 10 10 M Q , where the data are > 80% 
complete. Incompleteness affects mostly red galaxies, 
hence the plot shows > 80% of the SF galaxies. Their 
current SFR would generate their observed M t within td 
of 6 x 10 s to 6 x 10 9 , yr, mostly < 3 x 10 9 yr. These 
td are smaller than half of the age of the Universe at 
0.75 < z < 1. If we assume a high z/(~ 5), 85% of the 
galaxies have a SFR above the past average (b > 1), and 
57% a starburst (b > 2, Kennicutt et al 2005) 10 . At face 
value, this is contradictory; it should not be possible for 
a majority of galaxies to simultaneously undergo a star- 
burst, which by definition should occupy a short part of 
a duty cycle. 

However, these galaxies form the MS, and their SFR 
show a gradual decline on Gyr timescales since z ~ 1 
(Paper I), likely even since z = 1.4, not an enhancement 
of SFR in the z range shown in Fig. 2. These galax- 
ies therefore do not seem to simultaneously undergo a 
brief (< 1 Gyr) episode of elevated SFR on top of lower- 
level SFHs; instead, their high SSFR represent the early, 
strongly star-forming phase of a SFH that smoothly de- 
clines on Gyr timescales to z = 0. Their high SSFR 
(short td) imply (i) that this gradually declining epoch 
must form the bulk of their M*, and (ii) that the galax- 
ies must be observed < ltd after the onset of this epoch, 
suggesting Zf < 2 for > 60% of these galaxies; otherwise 
the produced M± would be higher than observed. 

These lines of evidence suggest that the observed high 
SSFR of many galaxies are not due to a periodic or irreg- 
ular burst, late in their SFH. Instead, many such galaxies 
seem to be observed shortly after their "initial burst" 
phase, the early stage in their predominantly smooth 
SFH which forms most of their M*. Moreover, the av- 
erage SSFR increases smoothly to lower masses, at all 
z. This points to a smooth dependence of the average 

10 Even if all our diagnostics overestimated SFR by a factor of 
2, the fraction of galaxies with b > 1(2) would still be 57(21)%. 



Zf on galaxy mass. Based on this evidence, we propose 
a scenario of "staged galaxy formation", where the aver- 
age onset of the major SF (zf) decreases smoothly with 
galaxy mass. This scenario achieves high SSFR without 
requiring that a large fraction of galaxies at any epoch 
are elevated in SFR (b > 1) or starbursting (b > 2). Al- 
lowing lower z / for a fraction of less massive galaxies is 
the only possibility to avoid this contradiction between 
burst fraction and duty cycle. The staged r models we 
use to approximate these SFHs parametrize both the de- 
cline timescales (r), and the onset (zf), of the main SF 
episode as a function of mass. 

The range of td in Fig. 2 shows that the staged scenario 
only requires a fraction of less massive galaxies to form 
later: the range of Zf must reach to lower z for less mas- 
sive galaxies. In addition, the model does not exclude 
some low level SF prior to the onset of the major SF 
episode, effectively the epoch of assembly. This allows 
it to be consistent with the presence of old (~ 10 Gyr) 
stars in many low mass local galaxies (e.g. Grebel 2004). 

A relation between the galaxy mass and the onset time 
of the dominant SF episode is observationally and theo- 
retically supported: see Heavens et al. (2004); Iwata et 
al. (2007); Thomas et al. (2005) for early-type galaxies; 
Feulner et al. (2005), Reddy et al. (2006) and Erb et al. 
(2006) report a systematic decrease of stellar age with 
M* up to z = 3 — 5. CDM structure formation provides 
a framework for a systematic relation between the domi- 
nant SF epoch and present-day galaxy mass (Neistein et 
al. 2006). Finally, insofar as "downsizing" means that 
a characteristic epoch of high SSFR occurs early in high 
mass galaxies, while at low z only low mass galaxies ex- 
hibit high SSFR (Cowie et al 1996), a delayed onset of 
major SF in less massive galaxies is a natural part of this 
process. 
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